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C^ . in the reduced density matrix elements for a coupled system involving a pair of 

kicked rotors as also one made up of a pair of Harper Hamiltonians. These dynamical 

\^ ' fluctuations are found to constitute a reliable indicator of the degree of chaos in the 

Cp ■ quantum dynamics, and are related to stationary features like the eigenvalue and 

eigenvector distributions of the system under consideration. A brief comparison is 

made with the evolution of the reduced distribution function in the classical phase 

p ■ space. 
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1 Introduction 



Quantum chaos deals, broadly speaking, witli tlie study of quantum systems 
liaving chaotic classical counterparts. The problem of distinguishing between 
the quantum behaviour of systems with chaotic and regular classical dynam- 
ics has been looked into for quite a long time and static features such as level 
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repulsion and eigenvector distribution have been clearly identified as signa- 
tures of quantum chaos. However, while classical chaos essentially manifests 
itself in the dynamics of the system under consideration, analogous dynamical 
signatures of quantum chaos are not easily identified (see, e.g., (|l|l3|3|)). This 
is essentially because of the fact that the quantum time evolution of a closed 
and bounded system must exhibit quasiperiodicity due to discreteness of the 
underlying energy spectrum. The same quasiperiodicity can be seen in aperi- 
odically driven system due to the phenomenon of dynamical localisation (j^lsl). 

Consequently, the search for dynamical features of quantum chaos has pro- 
ceeded along different lines. In an early paper, Peres (J6|) pointed out that the 
quantum evolution of classically chaotic systems may show distinctive insta- 
bilities under external perturbations. A number of recent studies (|7|; Ig la llOl ) 
suggest that such instabilities indeed constitute a basic characteristic of such 
quantum systems. Characteristic signatures in the time evolution of such sys- 
tems are also to be found when these are coupled to other systems ( llUiallSl ): 



the instabilities manifest themselves in the time series of the reduced density 
matrix (RDM) of the system under consideration. In all these studies, the 
distinctive quantum signatures are sought in the response of the system under 
consideration to external influence. 

An alternative approach, to be pursued in the present paper, is to look within 
the system itself and to explore distinctive features of time evolution of its 
subsystems (an early work with a similar approach is to be found in (jl^). 
Though an isolated closed system does not show any random features in time 
evolution, its subsystems do. Recent investigations reveal that the reduced 
entropy (or some entropy-like object measuring entanglement in subsystems) 
obtained from the RDM of a subsystem does show up distinctive features for 
systems possessing the stationary signatures of quantum chaos, as compared 
to the regular ones (|l5l : lla llTt llSl : |l9|). Thus, this entropy (and the analo- 
gous entanglement measure) may be looked upon as an indicator of quantum 
chaos. We shall present evidence in this paper that the decay of autocorrela- 
tions in the temporal fluctuations of the RDM is strongly related with the lack 
of correlations among the Hamiltonian matrix elements of the system itself, 
and will relate this finding to the stationary features of quantum chaos. We 
will also present evidence indicating that the correlation among the elements 
of the Hamiltonian matrix falls off as the corresponding classical motion be- 
comes more and more chaotic (in this context, see (|20|)). It may be worthwhile 
to mention that such dynamical manifestations of quantum chaos, especially 
those relating to the entanglement between subsystems, r aay be par ticularly 
relevant in the context of quantum information processing (|2ll: 1221: l23l : 1241 \2l 



This paper is organised as follows. First, we outline (Sec. 2) our basic approach, 
following the papers (|l5l : llQI ). In Sec. 3 we follow (|l9r i in briefly presenting 
results for spin systems with two types of system Hamiltonians, one being 



described by a random matrix, and the other being a regular Harper's system 
(for a related model see (|26l)). In the next section (Sec. 4) we present results 
for subsystems of more realistic interacting systems. The following sections 
(Sec. 5 and 6) deal with a few other aspects of the problem, and concluding 
remarks (Sec. 7). 



2 Basic Approach 



The basic approach has already been discussed in an earlier paper (il9f ) (see 
also (15)). A pure state of a closed system continues to evolve unitarily as a 
pure state. Thus the von Neumann entropy defined in Eq.(2a) below (or the 
linear entropy defined in Eq.(2b)) remains constant irrespective of the nature 
of the Hamiltonian. However, the reduced density matrix of a subsystem i.e. 
density matrix traced out over states of the complementary subsystem, does 
not evolve in a unitary or reversible way. This explains why the entropy of a 
system coupled to a heat bath increases with time and may show signatures 
of chaos (|l3t). As seen in fjlSl). the semi-classical entropy Sd = — J p\n p dp dq 
(where p is the distribution function for an ensemble) calculated from the 
reduced Husimi distribution function for a subsystem also shows the rapid 
loss of correlation for a chaotic system (see below). 

This idea serves as a clue for investigating the dynamical features of quantum 
chaos for an isolated system. We proceed as follows. 

If a system S be a composition of two subsystems (A and B) , the density 
matrix of S may be traced out over the states of B, giving the reduced density 
matrix (RDM) of A ( say pRj.The system S evolves by the system Hamiltonian 
H as 

p{t) = exp{-iHt/h)p{0) exp{iHt/h). (1) 

Now, if pit) is traced out over B for each point of time to get the RDM of A 
(p/j(t)) at these successive time instants, then one can study the fluctuations 
in the matrix elements of pR{t) from the resulting time series. The schematic 
diagram for the time evolution of the RDM is given in Fig.l. 

One could study the fluctuations of the individual elements of the RDM by 
themselves. However, for the purpose of presentation it is best to confine our- 
selves to the following suitable variables viz. the von Neumann entropy(S'yAr) 
and linear entropy(5'L) defined as 



SvN = -Tr{pR\npR) (2a) 

SL = l-Tr{pl). (2b) 



P(t=T) 




Reduction 



Pj^(t=0) Evolution Pj^t=T) 

Fig. 1. Evolution of density matrix (p) and its reduced component (p/j). 

Both of them are indicative of the RDM elements fluctuations as a whole. 
Among the two, we choose the linear entropy for the presentation of most of 
our numerical results. Both these quantities are indicators of the degree of 
entanglement between the states of A and B. 



3 Two Prototype Systems 



We first consider two spin systems, once again by way of brief recapitulation of 
(15|;|l9|). As a simple prototype for quantum chaos we consider an A^ x A^ real 
symmetric Hamiltonian matrix Hf. with randomly distributed matrix elements 
for which the stationary features predicted by RMT are conformed to, and then 
we compare this with a second Hamiltonian matrix Hr (of the same dimension) 
that is equivalent to the Harper Hamiltonian on a torus (Eq.3). In terms of 
the co-ordinate (g) and momentum (p) on the torus, the latter reads 



H, = 71 cos(27rp/P) + 72 cos(27rg/Q) , (3) 

the torus being of area PQ with A^ basic states so that h = PQ/2'kN ()29t l8fll). 

Each of the systems may be considered as an assembly of n-interacting spins 
such that 2" = A^(see f)26l )). With the Hamiltonian of the system S defined as 
above, we focus on a subsystem A made up of, say, p(< n) specified spins and 
look at the reduced density matrix given by (using Eq.l) 



pnit) = Tr('^-^)[exp(^)p(0)exp(^)], 



h 



h 



(4) 



where H = H^ or H.,. as the case may be, and where Tr^^ ^' indicates partial 
trace with respect to states of the remaining (n — p) number of spins. We then 
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Fig. 2. Autocorrelation of the Hamiltonian matrix as defined in Eq.5 for (a) He 
(rapid initial fall is displayed in inset) and (b) Hr (71 = 72 = 2.0); (exponential fits 
are shown with dashed lines). 



(A), 



calculate S\^ (t) 
from Eq.(??). 



l-Tr{pn')) 



A^) 



or, alternatively, S\,j(^{t) = -Tr(pB In p/j)) 



While ref. ()19f ) presents the numerical results and is to be referred to for de- 
tails, figures 2 and 3 indicate graphically the distinctive characteristics of He 
and Hr for n = 7, and p = 4. Of these. Fig. 2 (see caption) depicts the auto- 
correlation among matrix elements of H^ and H^. While the autocorrelation 
involves products like if^, iH^^^ ;_,_^, we look at the correlation in a direction 
'parallel' to the diagonal elements, namely 

AH{m)=J2H,,Hl^^^^^, (5) 

id 

as being the most characteristic for the systems under consideration. 

In order to estimate the rate at which the autocorrelation falls off in a direction 
parallel to the principal diagonal, we fit the initial falling portion with an 
exponential fall and obtain the correlation length Ic- One finds from the figure 
that, as expected, the autocorrelation falls off much faster for H^. as compared 
to Hr- Following ref (J20|) (see also (|27l)). we suggested in ref. (|l9r i that the 
degree of correlation among the Hamiltonian matrix elements can be looked 
upon as an indicator of the degree of regularity in the quantum dynamics. In 
the same work we also suggested that the dynamical manifestation of quantum 
chaos might be sought in the degree of randomness of the fluctuations of the 
RDM elements of subsystems of the system under consideration. The latter 
can be estimated from, say, the rate of fall of the autocorrelation of the linear 
entropy as a function of time, depicted in Fig.3 for H^. and Hr- The typical time 
for fall of correlation is denoted with the same symbol, i.e., Ic, and one finds 
that there is indeed a good correspondence between the divergences in the 
values of Ic for the Hamiltonian matrix and for the temporal fall of correlation 



of the RDM elements. 
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Fig. 3. Autocorrelation of the time series of Sl for spin systems when the system 
Hamihonian is (a) He and (b) H,.\ h = 0.592 for both the cases (initial falls are 
shown on a shorter time scale in insets and exponential fits are given by dotted 
lines) . 



To summarise, it appears worthwhile to look for a reliable indicator of quantum 
chaos in the temporal fluctuations of RDM elements of subsystems. These 
fluctuations are of a random nature for chaotic systems as compared to regular 
ones, as borne out by the 'toy Hamiltonians' H^ and Hr- The degree of this 
randomness reflects, in a way, the lack of correlation among the elements of 
the Hamiltonian matrix of the system under consideration. Other measures 
for estimating the degree of randomness in the temporal evolution of the RDM 
are to be found in (|l9l). 



4 Coupled Kicked Rotors and Coupled Harper Systems 



We now turn our attention to systems made up of more realistic subsystems 



in order to see if the observations made in ref. 
can be claimed to have general validity. 



and briefly outhned above 



First, we consider below a system made up of two kicked rotors coupled to- 
gether. The next subsection will deal with a system composed of two coupled 
Harper Hamiltonians. 



4-1 Coupled kicked rotors 



Kicked rotors are widely referred to in the context of quantum chaos (see e.g. 
d^ l28l ) ) . A single kicked rotor is described by the Hamiltonian 



P 



H=— + KcosqJ2Sit 



riT) 



(6) 



where K is the kick strength and r is the interval between two consecutive 
kicks. The classical dynamics is described conveniently in terms of the standard 
map, for which chaotic features are known to dominate for K > Kc ^ .97 (for 
r = l). 
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Fig. 4. Autocorrelation Ajj among evolution matrix elements defined in Eq.8 for 
the coupled kicked rotor system; kick strengths are (a) K = W and (b) K = 0.1 
(see text); r = 1.0, c^ = 2.0 for both cases; corresponding correlation lengths Ic are 
shown in the figure. 

We consider the composite system made up of two kicked rotors and described 
by the Hamiltonian 



^ rn ^ 



+K2Cosq2^6{t 



riT] 



Cr smgi smg2, 



(7) 



where Cr is the coupling strength. Denoting the RDM for one of the sys- 
tems by pr, we look at the autocorrelation of the time series for Svn = 
—Tr (pijln (pr)) for various values of Ki = K2 = K{saj). 

The single step evolution operator f/ is a complex two-dimensional unitary 



array, and the autocorrelation of the elements of this array (we consider a 
finite dimensional truncation) can be expressed in terms of 



Au(m) 



E 

k,l 



^k,l\ \^k-\-m,l+m\ 



where we have once again limited ourselves to correlation in a direction parallel 
to the principal diagonal (as in Eq. (5)) and have considered only the modulii 
of the matrix elements, ignoring the phases (the latter lead to similar results). 

Fig. 4(a,b) depict the variation of Au{m) with m for two values of the non- 
linearity parameter K (see caption), and for the coupling strength Cr = 2.0 
. The choice of parameters is such that for Fig. 4 (a) the classical phase space 
of each rotor as also the phase space of the composite system is dominated 
by chaotic orbits, while for Fig. 4(b) the corresponding phase spaces are pre- 
dominantly regular (note that for small values of both Ki and K2, the phase 
spaces are regular regardless of c^., since the Hamiltonian can then be effec- 
tively written as the sum of two uncoupled free rotors). One finds a close 
correspondence between the classical phase space and the correlation among 
the matrix elements of the single step evolution operator (the basis chosen is 
the direct product of the momentum states of the individual rotors). 
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Fig. 5. Autocorrelation of the time series for Svn (see text) taking one of the coupled 
kicked rotors as subsystem where parameters as in Fig. 4 (exponential fits in thick 
lines). 

Fig. 5(a), (b) depicts the temporal fluctuation of the RDM (obtained by taking 
partial trace over the rotor with variables (^2, P2)), expressed through the 
variation of autocorrelation of Svn as indicated above (parameter values are 
the same as in Fig. 4 (a), (b) respectively). One flnds a close correspondence be- 
tween the fall of temporal autocorrelation in Svn (similar results are obtained 
for Sl [not shown here]) obtained for a subsystem and the lack of correlation 
among matrix elements of the single step evolution operator U of the compos- 
ite system. This once again goes to show that the temporal fluctuations in the 
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Fig. 6. Temporal variation of autocorrelation of Sy-M when (a) K\ = {K2 = 10 
and (b) Ki = K2 = (c = 2.0 for both the cases; initial fluctuations in shorter time 
scale are in the insets, the fits are thick lines in the main windows ). 

RDM elements can indeed be looked upon as an indicator of quantum chaos. 

It is important to note that the randomness of the reduced density matrix 
fluctuations is a property of the (composite) system under consideration, and 
not of the subsystem whose RDM we are looking at. In other words, the 
nature of RDM fluctuations depends on whether or not the system as a whole 
is chaotic, and is random for any arbitrary choice of the subsystem. This is 
seen from Fig. 6 (a), (b), where, in (a) we take Ki = 0, i.e., the flrst rotor, for 
which the RDM is constructed, is regular, while K2 and Cr are given values 
such that the entire system is predominantly chaotic. We compare this with 
(b) where we take Ki = K2 = 0, so that the entire system, as well as both 
subsystems, is regular. The rapid loss of autocorrelation in (a) compared to the 
quasiperiodic fluctuations in (b) tells us that the RDM fluctuations are indeed 
characteristic of the system under consideration, and not of the subsystem one 
is looking at. 



4-2 Coupled Harper systems 



In continuation of the observation made in the last paragraph, we now consider 
a system consisting of a pair of Harper Hamiltonians coupled together. Here 
each individual subsystem is regular, while the nature of the classical phase 
space of the coupled system is determined by the strength of coupling. 

The phase space for each subsystem is periodic in both q and p with periods, 
say, Q and P respectively. In other words, the phase space can be taken as 
a toroid with dimensions Q, P. The quantum mechanics on a toroidal phase 
space imposes an additional condition PQ/2Tih = N, where N must be an 



integer (see ( 23: l3 



We take the coupled system as 



H = Hi + H2 + Ch sin gi sin ^2 



(9) 



where Hiand H2 are two Harper Hamiltonians (Eq. 3) with identical parameter 
values and c^ is the coupling constant. The form of the interaction has been 
chosen such that the phase space of the composite system is a AD torus while 
the quantisation condition is to be separately imposed for each set of q and p. 
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Fig. 7. Poincare section for the classical dynamics of the coupled Harper system 
when the coupling strength is (a) 10 and (b)O.l (Li = L2 = Ki = K2 = 2 for both 
the cases). 

Fig. 7 shows the classical phase space in a Poincare section along qi-pi where 
one finds that the degree of chaos in the composite system can be tuned 
through the coupling strength Ch- Fig. 8, on the other hand, gives the time 
series for autocorrelation of the RDM elements, expressed through the linear 




Fig. 8. Autocorrelation of the time series of Sl{= l — Tr{p'j^)) for one of the coupled 
Harper systems when the coupling strengths are as mentioned in Fig. 7; h = 0.628. 
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Fig. 9. Autocorrelation of the Hamiltonian matrix as defined in Eq.5 for the coupled 
Harper systems for coupling strengths as mentioned in caption of Fig. 7. 

entropy, for the same two values of the couphng strength. One finds a close 
correspondence between the two sets of data which show that a classically 
regular (resp. chaotic) phase space indeed results, in the quantum context, 
in a quasiperiodic (resp. random) variation in the RDM elements. Addition- 
ally, Fig. 9 depicts the correlation among the Hamiltonian matrix elements, 
once again for the same two values of the coupling strength, and indicates, 
as already observed, a close correspondence between the degree of correlation 
among these matrix elements and the nature of quantum dynamics as revealed 
through temporal fluctuations in RDM elements. 




Fig. 10. Power spectra of fluctuations in the time series of Tr{p'j^) for the coupled 
Harper systems with coupling strengths as mentioned in caption of Fig. 7. 

The degree of randomness in the RDM fluctuations is also revealed in the 
power spectrum of the time series for, say, Tr{p'j^), which is simply the Fourier 
transform of the corresponding temporal autocorrelation. A random fluctua- 
tion is indicated by a broad-band power spectrum while a more regular tem- 
poral variation corresponds to a power spectrum with a few localised peaks. 
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A convenient way to quantify this (|l5l : llOT l is to compute the participation 
ratio (PR) for the power spectrum. Fig. 10 depicts the power spectrum for 
the coupled Harper model with the same parameter values as in Fig. 8 and 9, 
and once again confirms our basic observation that the RDM fluctuations are 
a reliable indicator of the degree of quantum chaos charactering the system 
under consideration. 



RDM Fluctuations related to Eigenvalue and Eigenvector Statis- 
tics 



The evidence presented above seems to indicate that the degree of correla- 
tion among the Hamiltonian matrix elements of a classically chaotic system is 
characteristically small compared to that for a regular one (which is already 
implied in the random matrix theory (see (J20|;|27|)) and corresponds to the sta- 
tionary features related to eigenvalue and eigenvector distributions for chaotic 
systems), and this, in turn is correlated with the degree of randomness in the 
temporal fluctuations of RDM elements of subsystems. We indicate below how 
this correlation comes about. 

Meanwhile we take note of the fact that, in the classical situation, there exists 
a whole spectrum of the degree of randomness in the underlying dynamics, 
ranging from fully regular motion, through 'soft chaos', to motion charac- 
terised by 'hard chaos'. One can be more specific and quantify the degree of 
chaos through such indicators as the measure in phase space of the region 
occupied by chaotic trajectories, and the average decay time of autocorrela- 
tions. Our numerical results presented above indicate that a similar spectrum 
of the degree of chaos can also be identified in the quantum regime through 
measures related to the RDM fluctuations (see (|l9l ) for evidence based on a 
model Hamiltonian with a parameter effecting a convenient tuning for chaos), 
as also through ones related to stationary features such as the nearest neigh- 
bour level statistics (e.g., the Brody distribution that interpolates between 
the Wigner and Poisson distributions (jit 121: 1311)). Our analysis below tends to 
confirm these observations. 

The reduced density matrix elements are given by 



iPR)mn - zJ P"^'."' 
I 

where p is written in the usual block matrix form, and pmk,ni is the matrix 
element in the fc'th row and /'th column of the m-n'th block of the matrix. In 
all our computations we have started with an initial composite density matrix 
which is a direct product of density matrices for the subsystems chosen : 
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p = Pa® Pb- 



Now, according to Eq.4, 



ml,nl 
I 



ipR)m,n = J2p 



l{Ea - Eb)t' 



X 



h 

u,u \ / 

J2{m,l\ Ea) {Ea\ p{t = 0) IE,) {E, |n, I) =J2 exp (-icOabt) 0l, 



a,b 



where \Ea) and \Ei,) are energy eigenstates of the system under consideration, 
with energy values Ea and Eb respectively. In the last line, Uab is {Ea — Eb)/h 
and 

J2 {Ea\ Pit = 0) \Eb) i{Eb\ n, I) (m, l\ Ea)), 



m,n 

I 



the last factor being nothing but the overlap between the m'th block of \Ea) 
and n'th block of |-Eb). 

Thus, 



'm,n 



2 



Tr {pI) = E {(Pr) 

m,n 

= E fE<ntfl)x exp{-t{uJab+UJa'b')t). (10) 

a,b,a'b' \m,n ) 



This shows that Tr{p^p,) and hence the linear entropy S'^must have Fourier 
components with frequencies depending upon combinations of available energy 
intervals [Ea — Eh) (quasi-energies in the case of time periodic systems such 
as the kicked rotor). But the amplitudes of these Fourier components depend 
upon the overlap of energy eigenstates ( quasi-energy states for kicked systems 
) over all available blocks for the subsystems. Therefore, for Hamiltonians 
with a preponderance of large energy intervals and with eigenstates extended 
over all the basic states, the power-spectra of Sl will be broadband ones 
with considerably large Fourier components for a wide range of frequencies. 
In the opposite case of Hamiltonians with a narrow range of available energy 
intervals, there will be only a relatively small number of components spread 
over a narrow frequency range. Evidently, this distinction will also show up 
in the temporal autocorrelations of the RDM elements, since the latter are 
simply the Fourier transforms of the corresponding power spectra. 
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Available Energy Intervals 

It is known from the random matrix theory that the nearest neighbour level 
spacing distribution (NNLSD) of energy spectra for a real symmetric random 
Hamiltonian matrix is of the Wigner type, while it is of the Poisson type for 
a regular Hamiltonian matrix. 
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Fig. 11. Distribution of all possible energy intervals (divided by h ) for the Hamil- 
tonian (a) He and (b) Hr (parameters as in Fig. 2). 
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Fig. 12. Distribution of all possible energy intervals (divided by h ) for the coupled 
Harper Hamiltonian for coupling strengths as in Fig. 7. 

However, for our present purpose, the NNLSD is not of direct relevance, and 
of greater consequence is the distribution of energy intervals between all pairs 
{Ea, Eh) determining the distribution of frequencies Uab = {Ea — Eh)/h. Since 
the available energy intervals are nothing but additive combinations of near- 
est neighbour level spacings, one expects that in accordance with the Wigner 
statistics, the distribution of Uab will cover a wider range for a random Hamil- 
tonian as compared to a regular one. 
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This is corroborated in Fig.s 11 and 12 for the spin systems (Sec. 3) and the 
coupled Harper system (Sec. 4. 2) where one finds that the energy intervals 
for a system with a greater degree of randomness in the Hamiltonian matrix 
elements are indeed spread over a larger interval as compared to the corre- 
sponding intervals for a regular Hamiltonian of similar type. 

In summary, the Wigner distribution of nearest level energy spacings, with 
its attendant distribution of available energy intervals, explains the dynamic 
features of quantum chaos as revealed through the RDM fiuctuations. 



6 Hybrid Hamiltonians 



While the above can be looked at as a prima facie explanation of the char- 
acteristic features of RDM fiuctuations, we present below further evidence of 
the role of energy eigenvalue and eigenvector distributions in generating these 
fiuctuations, which may serve as a pointer towards a more detailed future ex- 
planation. For this, we consider a pair of contrived hybrid Hamiltonians with a 
view to explore separately the roles of eigenvalue and eigenvector distributions. 
We note that an arbitrary Hamiltonian H may be written as, 

H = V-^EV, (11) 

where V is the eigenvector matrix and E is a diagonal matrix with the eigenen- 
ergies as its diagonal elements. 

We first consider the spin systems discussed in Sec. 3, as prototype instances 
of regular and random systems. In accordance with Eq.ll, Hr and He may be 
written respectively as, 

Hr = VI ErVr 



and 



H, = V-^ErV,. 



Using these, we construct the hybrid Hamiltonians 

Hrc = V-^ErV, (12) 

and 

H,r = V-'E,Vr. (13) 

Evidently, Hrc has the eigenenergies of the regular Hamiltonian Hr and eigen- 
vectors of the chaotic Hamiltonian He- The case of Her is just the reverse, 
i.e., its energy eigenvalues are those of He and eigenvectors are those of Hr . 
We consider the time evolution for each of these two hybrid spin systems and 
compute the RDM for appropriately chosen subsystems, as in Sec. 3. 



15 



CD 

O 
CL 



30 

20 

10 





A 


I.Ot 
0.8 




B 


■ 


- 0.6 
1 0.4- 


1 




1 


0.2 


1 




L. 


n n 


, *HL.,_ 


40 80 120 


40 80 120 


Frequency 




Frequency 



2.0j 




0.5n 




-1 


1.5- 




C 


0.4 








D 




i 1.0 
o 

Q_ 






fe 0.3 
1 0.2- 






1, 1 




0.5 
n n 


^ 


If 


iJjJl 


0.1 
n n 


! 


1 i! ['} 
f\l\\\\\ 


1 , 1 







40 80 12 


c 


) 40 80 120 




Fn 


9q 


uenc 


y 




F 


-req 


uency 





Fig. 13. Power spetra of Sl{= 1 — Tr{p'j^)) for spin systems when the Hamiltonian 
is (A) Hr, (B) Hrc, (C) Her and (D) He (see text for explanation). 



Fig. 13 shows the power spectra (reported earher in ref . (l32l) and reproduced 
here for the sake of completeness) of the time series for Tr (p|.) with Hamiho- 
nians Hrc and Her, while those for Hr and He are also shown for comparison. 
One observes that the broadening is slightly less for Her as compared to He- 
This is due to the fact that though the available frequencies {ujab + ^a'b') are 
the same for the two cases (since these depend on the available energy intervals 
alone), the amplitudes Wmn^mn) (^^e Eq.lO) are small for most of the m, ra's 
due to the loacalised nature of the eigenvectors, especially for higher frequen- 
cies. Again, for Hre and Hr the available frequencies cover a small range as 
compared to those for He (see Sec. 5), resulting in narrower power spectra; still, 
the power spectrum is somewhat broader for Hrc as compared to Hr because 
of the extended nature of eigenvectors in Ve- 
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7 Cocluding Remarks: the Classical Picture 



While there exists a^emiclassical analog of the von Neumann entropy (e.g., 
the Wherl entropy (33) defined in terms of the Husimi distribution) one can- 
not, strictly speaking, define a classical entropy in terms of the phase space 
distribution function p{q,p). On the other hand, the concept of Kolmogoroff 
entropy applicable to a classically chaotic system does not have a direct ana- 
log for quantum systems. However, one can still seek an analogy between the 
time evolutions of the reduced density matrix for a quantum system and the 
reduced distribution function for a classical one (see below). 

The main point we propose to make in the present paper, as in a couple of ear- 
lier papers, is that while dynamical features of classical chaos of a system are 
to be found in the time evolution of the system itself, analogous dynamical 
features for a quantum system are to be sought not in the system dynam- 
ics, but in the dynamics of subsystems. Thus, autocorrelations in the classical 
phase space of a chaotic system die down in the long run, while corresponding 
autocorrelations in the Hilbert space vary quasiperiodically in the quantum 
case. On the other hand, autocorrelations for subsystems exhibit similar dy- 
namical features for the two situations, namely a decay in the long run. For 
a classically chaotic system, one infers this from the fact that that autocorre- 
lations die down in the full phase space itself. For a quantum system, though 
the autocorrelations for the system itself do not die down in the long run, 
those for the subsystems do, as seen from the evidence presented above. 

What happens in both the classical and quantum situations as one performs 
the reduction to a subsystem is, of course, a loss of reversibility. In the clas- 
sical case, while the evolution of the distribution function in the full phase 
space is reversible and volume preserving, it is at the same time characterised 
by the feature of mixing. The individual points belonging to an ensemble are 
transported to far away regions of the phase space, developing complex pat- 
terns made up of whorls, tendrils and lacunae (J3J), and this feature of mixing 
prevails on reduction to a subsystem, while at the same time there occurs a 
loss in reversibility. 

For a chaotic quantum system, on the other hand, it is the reduction itself 
that brings out the characteristic features of the system dynamics. The re- 
duced density matrix is significantly different as compared to a regular sys- 
tem, and shows characteristic fluctuations of a random nature. While we have 
presented above a preliminary explanation of this feature, a more detailed 
analysis remains pending. 

We conclude with two remarks on possible future work. One relates to the 
formulation of criteria for quantum chaos without reference to the classical 
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situation. Indeed, it seems desirable to have independent criteria f or q uantum 
chaos with a view to systems having no classical analogs (see (jlST l for an 
instance relating to the Baker's map), as also for characterising the dynamics 
of systems in the deep quantum regime where the semi-classical analysis bears 
no relevance (|35l). The present approach, based on RDM fluctuations appears 
worthwhile from this point of view. 

The other observation one may be interested in relates to possible implications 
in statistical mechanics. Consider, for instance, the problem of approach to 
the canonical distribution, where the system S under consideration interacts 
weakly with a heat bath H, making up a larger composite system C. One looks 
at the evolution of S, reducing from that of C. The conventional approach is 
to go over to the limit of infinite number of degrees of freedom and a quasi- 
continuous spectrum for the heat bath H in order to arrive at the equilibrium 
distribution (see e.g. (|35l)). An alternative approach would be to focus on 
the possible role of a chaotic interaction Hamiltonian between S and H in 
bringing about the equilibrium distribution for the reduced density matrix of 
S. A preliminary presentation based on this point of view is in preparation 
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